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Abstract The main purpose of this work is to present an interpretation method for 
injectivity test in a multilayer reservoir, based on new analytical solutions to the well 
pressure response. The developed formulation uses a radially composite reservoir 
approach and considers that the water front propagation may be approximated by a 
piston-like flow displacement. The reservoir is assumed to be laterally infinite and 
properties such as permeability and porosity may be different in each layer. The so- 
lutions were developed in the Laplace domain and then inverted to real domain using 
the Stehfest Algorithm. The proposed formulation was then validated by compari- 
son with a numerical flow simulator. Results showed a good agreement between the 
numerical simulator and the analytical model. Also, a sensitivity study was done by 
comparing the results of different scenarios varying oil viscosities and injection flow 
rate to assess how these properties affect the pressure and pressure derivative profiles. 
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Article Highlights 


e Pressure response in injectivity test 
e Multilayer and laterally infinite reservoir 
e Validation with a flow simulator 


1 Introduction 
The injection of water into oil reservoirs in order to extract relevant information re- 


garding the characteristics of the reservoir has become an interesting alternative to 
conventional well testing, given that it has many advantages. In addition to bringing 
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greater safety to the procedure the method helps to reduce CO? emissions to the en- 
vironment, since there is no need to gas/oil flaring which happens in a conventional 
well testing. However, it brings a new parameter to the problem since the system 
changes from a single-phase flow to a two-phase flow. 

Operationally, the conventional and injectivity tests are completely different, how- 
ever from the injectivity well testing we have enough information to perform the 
pressure behavior analysis. 

With this analysis it is possible to obtain different reservoir properties, according 
to what it is requested or necessary to the study of the reservoir (Neto et al., 2020a). 
The method consists of injecting water at a constant rate during a period of time. 
After that, the well is shut in and then it is possible to analyze the pressure curves to 
estimate properties of the reservoir. 

Solutions for injectivity test in the real field applied to a single layer and a mul- 
tilayer reservoir have already been developed (Bela et al., 2019). Moreover, there 
exist solutions in the Laplace domain to injectivity test in a single layer reservoir 
(Neto et al., 2020b). But, to the best of our knowledge, this is the first time a solu- 
tion for two-phase flow in a multilayer reservoir developed in the Laplace domain is 
presented. 

The pressure solution was obtained using a system of linear equations built in 
Laplace domain and brought to the real field using Stehfest’s Algorithm (Stehfest, 
1970) including the boundary and coupling conditions adopted to the scenario studied 
here. This work extends the single-layer solution developed by Neto et al. (2020b) to 
multilayer reservoirs and the behavior observed in the curves of pressure and pressure 
derivative was analyzed in different cases considering a variation of properties in the 
layers of the reservoir. Also a computational model was built to compare the results 
obtained. 

This work is divided into 5 sections. After the introduction the conceptual model 
is presented, followed by the implemented mathematical model. Then, there is the 
results section, containing the different cases analyzed, and lastly the conclusion. 


2 Conceptual Model 


The formulation of the problem considers a reservoir with two layers both with two 
regions, where the outer radius is the same in both layers and there is no formation 
crossflow. The well is assumed to be at the center of a circular zone in a cylindrical 
model (Nie et al., 2012). In this geometry there is a radial change in the properties as 
it passes from region | to region 2. 

A schematic illustration of the physical model is shown in Figures | and 2. 

The external boundary is considered infinite acting. Region | represents the flooded 
zone and region 2 correspond to the oil region, both considered homogeneous, in ad- 
dition the interface of the regions is flat and it has no thickness. The fluid mobility 
(k/) and compressibility are different on each side and each layer thicknesses are 
constant. Also there is no resistance to flow between the two reservoir regions. 
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Fig. 1: Schematic illustration of the side view of the reservoir model 





Fig. 2: Schematic illustration of the top view of the reservoir model 


In the actual reservoir, a variation of properties around the wellbore at the inter- 
face of the regions can be observed due to the rock-fluid interaction. However those 
alterations were not considered in the model presented here. 

For the mathematical formulation, the following hypotheses are considered: 


1. Two-phase flow without interaction between the fluids and a piston-like fluid dis- 
placement in which the water pushes the oil completely. 

. Isothermal flow with negligible gravitational forces and ignoring the impact of 
capillary forces. 

. Same initial pressure (pj) in all layers and small pressure gradients. 

. Reservoir is considered isotropic along each layer. 

. Single vertical injection well. 

. Constant flow rate during injection time. 
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7. Constant thickness (/) and porosity (@) on each layer. 
8. No flow considered in 6 and z-directions. 
9. The well is fully perforated along its extension. 


3 Mathematical Formulation 


Considering the reservoir model described in Section 2 a system of linear equations 
in Laplace domain is built to represent the differential equations of the water and oil 
phases, as well as the boundaries and coupling conditions. 

The initial condition (IC) reflects the pressure distribution at the initial time along 
the reservoir. The internal boundary (IBC) indicates how the water is injected dur- 
ing the test, which, in this case, will be a constant flow rate. The external boundary 
condition (EBC) is relative to the flow behavior at the reservoir extreme limit. As 
it is considered a laterally infinite reservoir, that is, the reservoir boundaries, repre- 
sented by the EBC, will not affect the well pressure response during the test time. 
The mathematical formulation was developed based on a consistent unit system. 

Considering a piston-like fluid displacement, the fluid mobility (Ai), hydraulic 
diffusivity (7;) and the derived fractional flow ({f/,) are defined by: 


n k jkr ji 
yh ard IE 1 
Gita: (1) 
k jkr ji 
a= (2) 
Hji 
1 
/ 
——— 3 
a eal ait i) 


The layer properties will be indicated by indices ji, where j represents the layer 
and i the region. Therefore, k; is the layer permeability, k;;; is the relative permeability 
of the fluid. 

Using the Buckley-Leverett formulation (Buckley and Leverett, 1942), the water 
front radius (rf) is defined by: 








in i, (t) 
rr j(t) = |r? + oi is (4) 


where vjnj, is the water volume injected into each layer. 


t 
vin C) =f ading(t) a 5) 
Thus, the following equations model the problem: 
Region |: 
PDE: 





0 (eae) = 1 dpj(nt) 6) 


Or or nj at 
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IC: 
P(t) = pi (7) 
IBC: Seat 
(ryt 
qj = —2mAjthj (« ae ~ ) (8) 
Region 2: 
PDE: 3 eer ee 
Pj2\nt - PjAnt 
ar (: ar ) ~ hp ot ©) 
IC: 
Pp2(nt) = pi (10) 
EBC: 
lim pjo(7,1) =0 (11) 


The relationship between the regions and layers are expressed by the coupling 
conditions. The pressure and flow-rate equalities in each region at the interface be- 
tween them, that is, at the water front radius (rf), represent the coupling conditions 
by region (CCR) (Nie et al., 2011). This condition is defined as: 

CCR: 

Apj(trjt) =App(rr;.t) 
(12) 
qj(rrj,t) = 4j2(rFj,t) 


Applying Darcy’s Law to compute the flow-rate in each region: 
ny IP j\ 4 IP; 
2nd 1h; (+ AG ) = 20h johj (: > 


As one can see, both, water and oil rates are considered the same, however it is 
only an approximation that can be applied in a piston-like problem. 


Therefore, 5 
Op; A; Op; 
(: a) = Ape (: 52 
or THIF] Aji or TIF] 


Equations (13) and (14) refer only to the flow-rate CCR (12). 

Since the considered reservoir model assumes that layers have no connection to 
each other except by the wellbore, the coupling condition by layers (CCL) is given 
by the pressure equality and the principle of mass conservation at this same point as 


shown below: 
CCL: 








(13) 








THIF j THIF] 


(14) 














Apj(tw,t) = APj+i1(%w,t) 
n 





1 
dinj = X ajiltw.t) ( 5) 
j= 
Bringing the equations to Laplace domain: 
Region |: 
PDE: ore : 
Pji 2 
a i —pj(rt=0 16 
ry (Pp) = qowPn lew) —palet =0) (16) 


6 Ana Luiza Mastbaum et al. 





Applying the Laplace transform in equations (6) to (11) and applying the IC 
(given by equations (7) and (10)), the following Ordinary Differential Equation for 
Region | is obtained: 











ODE: 
re) Py(rnu) 1 Opi (r,u) 1 
= Dj — Pi 17 
a Pe. ae na (upj(%4) — pi) (17) 
IBC: : 
Apart —2A jh; 1 
(2) ’) = oN (18) 
Or J\romy, jl 
Region 2: 
ODE: , 
O-pp(rnu)  1dpp(ru) 1 
= Dj — Pi 19 
ae a ee As (up j2(%,4) — pi) (19) 
EBC: 
lim p ja(r,u) =0 (20) 
The pressure change in Laplace domain may be defined as Ap(r,u) = p; — p(r,u) 
and oe = ane. Thus, applying this definition for pressure change in equations (16) 
to (20), the following equations are obtained: 
Region |: 
ODE: 


PAP 1A AP;  u 


























Apj, =0 21 
or ror nj ee 
IBC: 3 
Ap; —qj 1 
(« Pj ) = Gil er (22) 
or r=ry 2nAj1h; Uu 
Region 2: 
ODE: 5 : 
fe) Apj2 lo Apj2 u 
Apj2=0 23 
or ror nj2 i a 
EBC: 
lim ABj(7,u) =0 (24) 
Besides that, the coupling conditions in the Laplace domain are shown as: 
CCR: 
Apji(rrj,t) = AP ja(rrj,t) 
AAPA _ Ap ( AAP? (25) 
Or aa at 
TSIFj TIF j 
CCL: 
APj,1 (Tw,t) = APj+1,1 (Mw, t) 
inj a e (Se ) (26) 
: or 


nN» 
2mu X Aiki j=l r=ry 
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The pressure solution for any layer j and region i in Laplace domain is: 


AG = ALR (« =| Bele (+ =| (27) 
V ng; V nei 


where A ;; and B;; are coefficients to be determined. It is also known by the Bessel 
functions properties (Abramowitz and Stegun, 1964) that: 


lim Ko(r) =0 (28) 
roo 
lim Io(r) = 0 (29) 
r—yoo 


Thus, to satisfy the external boundary condition of Region 2 given by equation 
(24), Bj2 must be zero. Consequently, the pressure solution for i = 2 is: 


Apne (- =) (30) 
1j2 


By the IBC (22) and the flow-rate CCL (26): 
n ( | 
3 Ger 
= or 
j 


Thus, deriving the equation (27) with respect to r = ry, may be combined with 
equation (31): 


n Uu Uu 1 inj 1 
i: (anki (rm / =) — Bah (" )) = = x = (32) 
j=l Nj Nj Tw na on Y Ajuh)) Uu 

j= 


Replacing the pressure solutions for Region 1 (27) in the corresponding CCL 
(26): 


u Uu Uu u 
A\Ko & / =) + Bylo (rm / =) = A21Ko (> (=) + Boilo ( =) 
jl Njl Nl jl 
(33) 


On the other hand, from equations (25) and (27), the following expression for the 
flow-rate CCR is obtained: 


u u ae Nil ( =) 
Aki ( re;,/— )—Byh ( rr;,/ = fA pK, ( rej,/— 34 
at (re 7) at 1G a) Aj V 172 ee etd 1j2 & 


The last condition to be applied to this model is the pressure CCR (25). Therefore, 
combining the pressure solutions (27) and (30): 


Aj1 Ko (v«) ~) + Bylo (r+) =) = Aj2Ko (v«) =~) (35) 
V nA Vina j2 


inj 1 
See rae a) 
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The equations of hydraulic diffusivity, boundary conditions and coupling condi- 
tions of each region and layer provide a linear system that models this problem in the 
Laplace domain. The following equality represents this system: 


1 


qinj 





rite 2 u 
All wal yy 20 Dy (Ajihj) 


By 
Ai2 
A21 
By, 
Ax 


=1 


(36) 


oooco co. 


Where M is a linear system in matrix form composed of equations (32), (33), 
(34) and (35). Equation (32) represents the coefficients of the well pressure equation. 
Equation (33) expresses the coupling condition between layers and equations (34) and 
(35) express, respectively, the flow-rate and pressure coupling conditions between 


regions of both layers. 


For a better understanding, the matrix M was divided into two parts, as shown 
below. The first part consists of its first three columns and the second one by its last 


three. 














Ko |rriy/7,,) 40 (re1/4) 
0 0 
0 0 





Ky (rs x) —i (i [-) 
Ko CG 7 ) Io (ah [-) 


Ka (r/c) A (reves) ~ 
(rr (ae) 





dot / Twi fe 
wl Nol Ky (re Nol ) (37) 
—Ko (re / x) 


0 





0 


Thus, the last three columns of matrix M are represented by: 
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0 0 0 
—Kp (rs Vis) —I la j=) 0 
0 0 0 
0 0 0 i 


A 
Ki (reay/aia) Wi (reay/ a) 2 REA (122 a) 
Ko (rr2/5) Ip (rr2\/5) —Ko (rr2,/5) 


Therefore, once the coefficients A;; and B,; were determined, it is concluded that, 
in Laplace domain, the well pressure solution is calculated by: 


APwr =A11Ko | rw ane Bile tyne (39) 
Tw Nw1 


The well pressure profile given by equation (39) is then converted to real domain 
through Stehfest Algorithm (Stehfest, 1970). 











4 Results 


The accuracy of the analytical model detailed in Section 2 was determined by com- 
parison with a commercial flow simulator (von Hohendorff Filho, 2004). The curves 
of pressure (AP) and pressure derivative with respect to the natural logarithm of time 
(AP’) were plotted on log-log graphs functions of time. After the proposed solution 
was validated, a sensitivity study was performed, evaluating the influence of oil vis- 
cosity and injection flow-rate in the computed pressure response. 


4.1 Model Validation 


The endpoint mobility ratio is a ratio between the mobility of the fluids under discus- 
sion. It is a critical factor to influence flooding efficiency of the water in the reservoir. 
When the endpoint mobility ratio is greater than one, it means the water is more mo- 
bile than oil. Otherwise, it means the water is less mobile than oil, leading to a better 
sweep efficiency. Thus the mobility ratio is defined by: 


wee (40) 
Ao Lw 
The results are divided in two cases, considering different mobility ratios between 


the fluids. The reservoir properties for both cases as well as the injection flow-rate are 
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Table 1: Input data for both cases 





Parameter Value Units 





inj 800 m3/day 
ky 125 mD 

ky 250 mD 

hy 15 m 

hy 15 m 

On 0.25 fraction 
o2 0.25 fraction 





Table 2: Fluid properties in each case 





Case Uy (cP) by (cP) M 





1 5.0 10 
2 0.2 ae 0.4 





shown in Table 1. Meanwhile, the fluid viscosities for Cases 1 and 2 can be seen in 
Table 2. Also, it was considered a 96-hour flow period on both cases. 

In Case 1, a 5 cP oil viscosity was tested, leading to a high mobility ratio and a 
better movement of the water, although not a great sweep efficiency in moving the oil 
along the radius of the reservoir. In Case 2, a lower oil viscosity was tested, leading to 
a mobility ratio lower than in Case 1 and consequently a better scenario to the sweep 
efficiency. 

Figures 3 and 4 show the pressure and derivatives curves considering both the an- 
alytical solution and numerical data. A slight deviation between the curves at the short 
time can be observed due to the simulation grid refinement causing this line spacing 
can be observed. Besides that, in both cases it is possible to notice the proximity of 
the results showing a good agreement between the models. 

In addition, this refinement problem appears to be more sensitive in Case 1, even 
for a long time. It can be explained by the water front behavior, i.e., piston-like and 
due to its mobility ratio, since it is higher than one and greater than Case 2, leading 
to a more unstable displacement. 

It is also observed that the derivative starts at a certain level and then stabilizes 
at another one, which is the same for both cases, even though the oil viscosity is 
different. This can be explained by the fact that, in an injectivity test, the first level of 
the pressure derivative curve represents the oil mobility, while the second represents 
the water mobility. Thus, since the properties of the injected water are the same for 
Cases | and 2, the second derivative level is expected to be the same in both graphs. 


4.2 Sensitivity Analysis 
In order to analyze the proposed solution with different scenarios, the pressure and 


pressure derivative over time were plotted, in five different oil viscosities (maintain- 
ing the injection rate at 1400 m?/day) as shown in Figure 5. 
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Fig. 3: Pressure and Pressure Derivative data for Case | 
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Fig. 4: Pressure and Pressure Derivative data for Case 2 


The pressure derivatives for different oil viscosities in Figure 5 converge to the 
same level, as well as Figures 3 and 4. The reason behind this behavior is the same 
as mentioned in Section 4.1. As the oil viscosities are different for each curve and 
the first level of the pressure derivative curves represents the oil mobility, their initial 
level are not the same. However, as the properties of the injected water remain the 
same, the second level of all curves converge to a single threshold. 
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Fig. 5: Pressure and Pressure Derivative for different 011 viscosities 
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Fig. 6: Pressure and Pressure Derivative for different injection flow-rates 


It is also possible to notice that the pressure curves rise as the oil viscosity in- 


creases. This is due to the fact that the higher the oil viscosity, the more distant it is 
from the water viscosity and therefore it departs from the 0.5 cP level. When oil and 
water viscosities are the same, pressure derivative remains constant during the entire 
test, since there is no difference between fluid mobilities. 


Figure 6 displays the variation of injection rate maintaining the same oil viscosity 


as Case |. Figure 7, on the other hand, enables easy visualization of the difference 
between the behavior of each scenario, once it shows the normalized pressure and 
pressure derivative data. The normalized pressure is given by each pressure divided 
by its injection rate versus time. Thus, pressure derivative is also normalized by the 
injection flow-rate. 
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Fig. 7: Normalized Pressure and Pressure Derivative for different injection flow-rates 


The decrease of the derivative pressure curves indicates the transition region be- 
tween the oil properties region and the water properties region. Thus, looking at those 
same curves in Figure 7, it is noticed that as the injection flow-rate decreases it takes 
longer to reach the water level. The opposite occurs as the injection flow-rate in- 
creases. This is due to the fact that, with a lower injection rate, it takes longer for 
the water bank around the wellbore to become relevant, which means that pressure 
change will reflect water properties. However, considering the normalized pressure 
derivative data, all curves achieve the same level once the transition is over. 


5 Conclusion 


Based on the solution for pressure change during injectivity tests in single-layer reser- 
voir, it was possible to develop a formulation to a multilayer reservoir. The proposed 
multilayer solution uses coupling conditions between layers and between regions to 
build the linear system that must be solved to determine the pressure change. 

The comparison between the numerical model and the analytical model developed 
in this work shows the proximity of the results and accuracy of the solution, once both 
tested cases obtained a good agreement with the numerical simulator. 

The sensitivity study was done comparing the results considering different oil 
viscosities and injection flow-rates. From this sensitivity study, it was possible to 
notice how these parameters influence the pressure response and pressure derivative. 

Thereby, this work presents a new and accurate solution for pressure response 
during a water injection well testing in a way that is possible to evaluate the wellbore 
pressure response. 

Although the work was based on a reservoir with only two layers, the presented 
solution may be extended to an arbitrary number of layers. This interpretation was 
possible considering an infinite reservoir with the boundary conditions not affect- 
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ing the pressure behavior. However the model can be adapted to different boundary 
conditions, for that the definition made on EBC (11) must be modified. 
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Nomenclature 

Pp pressure in Laplace domain 

N ji hydraulic diffusivity of the fluid in region i and layer j 
Ai fluid mobility of the fluid f, mD/cP 

M mobility ratio, fraction 

i viscosity of the fluid in region i, cP 

@ porosity, fraction 

Te fractional flow derivative 

hj layer thickness, m 

I, modified Bessel function of first kind and order v 

kj layer permeability, mD 

K, modified Bessel function of second kind and order v 
ky relative permeability of the fluid in region i, mD 


Di initial pressure, kgf/cm” 

Dj layer pressure, kgf / cm 

Pw bottom-hole pressure, kg f/ cm 

di flow-rate of the fluid in region i, m? /day 
inj injection rate, m> /day 
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Tw wellbore radius, m 

IF j water front radius in layer j, m 

Sox residual oil saturation, fraction 

Si initial water saturation, fraction 

t time, Ars 

Vinjj water volume injected into each layer, nm 
c compressibility, (kgf/cm?) ~! 
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